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As an example of an application of lattice QCD we describe a computation of four-quark 
operators in the nucleon. The results are interpreted in a diquark language. 

1. INTRODUCTION 

Lattice gauge theory offers the possibility to perform non-perturbative computations in 
a strongly coupled theory without any model assumptions. Lattice QCD, in particular, 
allows us to investigate low-energy properties of QCD, like, e.g., the hadron mass spec- 
trum, decay constants, form factors, and moments of structure functions. In this article, 
we shall discuss as an example the computation of four-quark operators in the nucleon 
recently performed by our group (the QCDSF collaboration). Since four-quark operators 
are related to higher-twist effects in nucleon structure functions, a subject of considerable 
theoretical and experimental importance, they are a particularly interesting object for 
lattice-QCD studies. For more technical details we refer to Ref. 111]. 



2. WHAT CAN WE COMPUTE ON THE LATTICE? 

The basic observables in lattice QCD are Euclidean n-point correlation functions. Since 
space-time has been discretised (with lattice spacing a) the path integral has become a 
high-dimensional ordinary integral over a discrete set of field variables, which can be eval- 
uated by Monte Carlo methods. As the (Grassmann valued) quark fields appear bilinearly 
in the action, they can be integrated out analytically leaving behind the determinant of 
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the lattice Dirac operator and products of quark propagators. In the quenched approxi- 
mation, which will be employed throughout the whole paper, this determinant is replaced 
by 1. The quenched approximation saves a lot of computer time, but it is hardly possible 
to estimate its accuracy. 

Let us briefly sketch the computation of hadronic matrix elements. First, one has 
to choose suitable interpolating fields for the particle to be studied. For a proton with 
momentum P one may take 

B^{t,P)^ e-^%,<(x)4(x)(C'75)/37<(x) (1) 

and the corresponding B, where i, j, . . . are colour indices and a, /?, . . . are Dirac indices. 

As the time extent T of our lattice tends to oo, the two-point correlation function 
becomes the vacuum expectation value of the corresponding Hilbert space operators with 
the Euclidean evolution operator e~^^ in between, i.e. we have, omitting Dirac indices 
and momenta for simplicity: 

(5(t)£(0)) ^=°° (0|fic-^*fi|0) . (2) 

If in addition the time t gets large, the ground state \p) of the proton will dominate the 
sum over intermediate states between B and 5, and the two-point function will decay 
exponentially with a decay rate given by the proton energy Ep-. 

{B{t)B{0)) (0|Be-^*B|0) * (0|B|p)e-^-*(p|B|0) + • • • (3) 

Of course, if the momentum vanishes, we have Ep = nip, the proton mass. 

Similarly we have for a three-point function with the operator O whose matrix elements 
we want to calculate: 

{B{t)O{T)B{0)) (0|Se-^(*-")Oe-^^S|0) 

= (0|E|p)e-^^(*-")(p|C|p)e-^->|S|0) + • • • 

= {0\B\p)e-^^\p\B\0){p\O\p) + --- (4) 
if i > T > 0. Hence the ratio 

should be independent of the times r and t, if all time differences are so large that excited 
states can be neglected. 

The proton three-point function for a two-quark operator contains quark-line connected 
as well as quark-linc disconnected pieces. In the quark-line connected contributions the 
operator is inserted in one of the quark lines of the nucleon propagator, while in the 
disconnected pieces the operator is attached to an additional closed quark line which 
communicates with the valence quarks in the proton only via gluon exchange. In the limit 
of exact isospin invariance as it is considered in this paper, the disconnected contributions 
cancel in the case of non-singlct two-quark operators. Moreover, one might try to argue 
that it would be more consistent with the quenched approximation to drop quark-line 
disconnected contributions generally. 
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3. WHAT DO WE WANT TO COMPUTE ON THE LATTICE? 

The operator product expansion (OPE) relates moments of structure functions to nu- 
cleon matrix elements of local operators. In the deep-inelastic limit, where the momentum 
transfer becomes large, one can express Nachtmann moments of, e.g., F2 in the fol- 
lowing form: 

ll dxx--'F,{x, Q^)U,,,,_ = c(f)Af (^) + ^^'H/i) + O (^) • (6) 

Here n = 2, 4, 6, . . ., the (reduced) matrix elements A^*^ of twist t and spin n are renor- 
malised at the scale fi, and c^^ = c^^Q"^ / fj^, g{jj,)) are the Wilson coefficients, calculated 
in perturbation theory. Note that we go back to Minkowski space in this section. 

Besides the leading twist-2 contribution we have included in (^) the twist-4 correc- 
tions. These are suppressed by and have been the subject of intensive investigations 
both experimentally and theoretically. In the following we shall describe what lattice 
computations can tell us about twist-4 effects. 

An important class of twist-4 operators consists of four-quark operators. In particular, 
the twist-4, spin-2 matrix element A^2^ is defined by the following expectation value in a 
proton state with 4-momentum P ({■■■} indicates symmetrisation) 

{P\Al^^y - traces|P) = 24^^(P^P^ - traces) (7) 
with the four-quark operator 

A% = ^PG^,^,t^^P^PG^,^,t^i;. (8) 

Here the quark field ip carries fiavour, colour, and Dirac indices, the are the usual 
generators of colour SU(3), and G is a matrix in fiavour space, the charge matrix 

G = diag(e„, ed) = diag(2/3, -1/3) . (9) 

(2) 

The corresponding leading-twist matrix element A2 is defined by 
(P|C|^^} - traces|P) = 2Af^(P/,P^ - traces) (10) 
in terms of the two-quark operator 

0^,= '-ijG^^^D,^. (11) 

This operator has dimension 4, while the four-quark operator (H) has dimension 6. Hence 
four-quark operators may mix with operators of lower dimension (two-quark operators). 
This is a rather general phenomenon for operators of higher twist, because for given spin, 
higher twist means higher dimension. Unfortunately, the mixing with operators of lower 
dimension is hard to deal with in perturbation theory. For the time being we do not 
attempt a non-perturbative calculation of the renormalisation and mixing coefficients. 
Instead we look for four-quark operators for which mixing with two-quark operators is 
prohibited by fiavour symmetry and apply perturbative (one-loop) renormalisation. 
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4. FOUR-QUARK OPERATORS IN THE PROTON 

Assuming SU(2) flavour symmetry (isospin symmetry) we find that two-quark operators 
can have isospin / = or 1, whereas four-quark operators can have / = 0, 1, or 2. Hence 
four-quark operators with 1 = 2 cannot mix with two-quark operators. Unfortunately, 
the expectation value of any 1 = 2 operator in the proton vanishes. So we have to enlarge 
the flavour symmetry from SU(2)f to SU(3)f, i.e. we assume three quarks of the same 
mass, and the flavour structure of the operator in the OPE is now 

= {cuUU + Cddd + esSs){euUU + Cddd + CgSs) . (12) 

Whereas two-quark operators transform under SU(3)f according to 3 ® 3 = 1 © 8, we 
have for four-quark operators: (3 (g) 3) (g) (3 (g) 3) = 2 ■ 1 ©4 ■ 8 © 10©T0© 27 . Four-quark 
operators with / = 0, 1, /a = 0, and hypercharge Y = from the 10, 10, 27 multiplets 
do not mix with two-quark operators and can be used in a proton expectation value, e.g. 
the 1 = 1 operator 

02j^ = —[{uu){uu) — {dd){dd) — {us){su) — {su){us) + {ds){sd) + {sd){ds) 

— {ss){uu) — {uu){ss) + {ss){dd) + {dd){ss)] (13) 

belongs to the 27 multiplet. Such a flavour component can be projected out by suitable 
linear combinations of the matrix elements of the octet baryons p, n, A, E, H, e.g. 

- 2(E°|0|E°) + = {p\Oif\p) . (14) 

Note that the flavour combination in four-quark operators which do not mix with two- 
quark operators (like £^2^^) is such that all quark-line disconnected contributions cancel. 
However, mixing with other four-quark operators (with different Dirac or colour struc- 
tures) remains and is taken into account using perturbatively calculated mixing coeffi- 
cients. 

After a chiral extrapolation (linear in the quark mass) we obtain in our special flavour 
channel 

1 dxF2(x,Q^)iyi=l,„ = -0.0005(5) ^ ^^'^ ^ + O(a^) . (15) 

The 0{al) contribution arises from the fact that we know only the tree level Wilson 
coefficient c^^ = g"^ [1 + 0{g^)) [0]. For reasonable values of the result is rather small 
using, e.g., as(Q^) ~ 0.3 for = 4GeV^. But it is interesting to note that the bag 
model estimates the prefactor in (|T5|) to be cx B/rrip ^ 0.0006, where B ^ (145 MeV)'' is 
the bag constant 0. 

5. DIQUARKS 

The four-quark operators can be rewritten to look like diquark densities. We have 
computed matrix elements of operators of the following form: 
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— {UaT'J:iCuJ){u^,C ^'j5^'Ub'){Sab'Sba' + Saa'hb') , (17) 

where a, b, a', and b' are colour indices. The Dirac structure of the operators is represented 
by the matrices F and T'. In (|TBp the diquark is in a 3 of colour and thus anti-symmetric 
in the colour indices. In ( p^Tf ) it is in a 6 and symmetric in the colour indices. Strictly 
speaking, we are again studying operators from the 27 representation of SU(3)f whose 
uuuu component is as given above. But at least within the quenched approximation 
the interpretation as valence diquark densities seems reasonable. In order to interpret our 
results we have combined the operators such that they correspond to diquarks of spin zero 
and spin one. For an operator Of^^, with two space-time indices (in Euclidean notation), 
e.g. for the Dirac structure F ® F' = 7^ (g) j^, we take the expectation value of O^i in 
a state with vanishing momentum to represent a spin-zero diquark and the expectation 
value of J2i=i C^ii to correspond to a spin-one diquark. 

In Fig. |l| we show our results for the diquark densities. They are plotted versus 1/k 
and linearly extrapolated to the chiral limit at 1/kct, where the variable 1/n determines 
the bare quark mass rrig = [l/n — !/«„.) /(2a). In physical units these bare quark masses 
are approximately 100, 190, and 250 MeV. 

The pattern of the results can tentatively be understood in a non-relativistic quark 
picture. When the diquark is in the 3 of SU(3)c it is anti-symmetric in the colour indices, 
and therefore the symmetric (in the spin indices) spin-one state is favoured over the anti- 
symmetric spin-zero state. On the other hand, when the diquark is in the (symmetric) 6 
of colour one might at first sight expect the anti-symmetric spin-zero state to dominate 
over the symmetric spin-one state. Although the spin-zero contribution is indeed less 
suppressed than in the 3 case it is not really dominating. This is probably related to 
the fact that a diquark in the 6 of SU(3)c must be accompanied by (at least) one gluon 
if it is to form a colour singlet with the remaining quark. The coupling to the gluon, 
mixing "large" and "small" components of the quark spinors, would invalidate the above 
arguments which worked reasonably well for diquarks in the 3 of colour. 
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Figure 1. Diquark densities with the Dirac structures (from top to bottom) F ® F' 
7//75 ® 7i/75 (3 of colour), cr^^ (g) a^^ (3 of colour), 7^ (g) 7^ (6 of colour). 



